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We study theoretically the fluctuation properties of acoustic phonons created in a semiconductor 
quantum dot after ultrafast optical excitation. An excitation with a single ultrafast pulse creates 
an exciton confined to the quantum dot, which is coupled to longitudinal acoustic phonons. This 
leads to the formation of a polaron in the quantum dot accompanied by the emission of a phonon 
wave packet. We show that the fluctuations of the lattice displacement associated with the wave 
packet after a single laser pulse excitation in resonance with the exciton transition are always larger 
than their respective vacuum values. Manipulating the exciton with a second pulse can result in a 
reduction of the fluctuations below their vacuum limit, which means that the phonons are squeezed. 
We show that the squeezing properties of the wave packet strongly depend on the relative phase 
and the time delay between the two laser pulses. 



I. INTRODUCTION 

The creation and manipulation of nonclassical quan- 
tum states of bosonic systems continue to attract large 
interest. Prominent examples of such nonclassical states 
are squeezed states. In a squeezed state the fluctuations 
of a given variable fall below their corresponding vac- 
uum value at the cost of increased fluctuations of its 
conjugate variable in order to satisfy Heisenberg's un- 
certainty relation. A well established field since many 
years is squeezed light, which can be generated in non- 
linear optical processes like parametric down conversion^ 
and which has applications in optical communication and 
measurements j 2 -^ 

Phonons are another type of bosons where nonclassi- 
cal states have become of growing interest in the past 
years. Squeezed phonons have been the subject of many 
experimental 4 8 and theoretical^"— studies. In most 
cases the considered phonons had a fixed frequency, ei- 
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ther because they were optical phonons 
Hove singularity appeared in the spectrum of acoustic 
phonons X An indication for squeezing is then the appear- 
ance of an oscillation with the double phonon frequency, 
however this alone is not yet an unambiguous proof for 
squeezing 

Phonons with a fixed frequency do not travel in the 
crystal because they have a vanishing group velocity. 
Therefore, squeezing produced in such kind of phonon 
systems remains in the region where it has been gener- 
ated. In contrast, in the case of photons squeezing can be 
transported at the speed of light from the place where it is 
generated to other places. Such a transport of squeezing 
could also occur for squeezed states of acoustic phonons 
which, due to their approximately linear dispersion rela- 
tion travel through the crystal at the speed of sound. 
An example for squeezing of acoustic phonons, which 
has been studied in the past, is the phonon analog of 
parametric down conversion, i.e., the decay of a coherent 



optical phonon into a pair of acoustic phonons Inter- 
estingly, however, it turned out, while indeed squeezing 
of the lattice displacement can be achieved in this case, 
the squeezing again remains where it has been generated 
and does not travel even though the generated phonons 
lead to a finite transport of energy away from the dot. 
The reason is that here the squeezing is strongly related 
to the quantum correlations between phonons with oppo- 
site wave vectors generated in the decay process. In this 
paper we will study a different process, the generation 
of acoustic phonons associated with the ultrafast optical 
excitation and manipulation of an exciton in a quantum 
dot (QD) structure. 

Quantum dots with their discrete level structure of 
electronic excitations are often referred to as artificial 
atoms. Due to their embedding in a surrounding semi- 
conductor matrix, however, they are much more strongly 
coupled to the environment than real atoms. Most im- 
portant for the dynamics of the QD exciton is in many 
cases the coupling to longitudinal acoustic (LA) phonons 
via the deformation potential interaction, leading to var- 
ious, often undesired phenomena like a phononic back- 
ground in absorption or luminescence spectr a 16 ! 17 or a 
damping of Rabi oscillations* 1 ^ — However, this coupling 
also gives rise to interesting spatiotemporal dynamics of 
the generated phonons, which typically consist of a lo- 
calized part remaining in the region of the QD and a 
traveling part leaving the QD. This latter part may be 
reflected at a surface and reenter the QD 21 or may travel 
to another QD influencing its optical properties^ Also 
bolometric measurements can be used to monitor the 
spatio-temporal dynamics of phonon wave packetS i 23 i 24 

In previous studies we have analyzed the possibility of 
generating phonon squeezing in the case of a QD cou- 
pled to longitudinal optical (LO) phonons i 11 ^ 2 We have 
shown that a single ultrafast excitation resonant to the 
exciton transition cannot create squeezing, while a reso- 
nant two pulse excitation can lead to squeezed phonons. 
For LO phonons with a single frequency it is possible to 
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transform the coupling to a coupling with a single ef- 
fective phonon mode. The dynamics can then be well 
illustrated by means of the Wigner function^ Such a 
reduction is not possible in the case of acoustic phonons 
with their continuum of phonon frequencies. It is the aim 
of this paper to analyze the fluctuation properties of LA 
phonons generated after ultrafast resonant excitation of 
the lowest QD exciton. In particular, we will show that 
indeed squeezed LA phonons may be created and that 
squeezed phonons can be emitted from the QD in the 
form of a phonon wave packet. 

The paper is organized as follows. In Sec. |TT] we in- 
troduce the theoretical model of the QD interacting with 
light pulses and LA phonons and we define the relevant 
variables for the study of phonon squeezing. Section ITTT1 
is then devoted to the results of our calculations, where 
we first consider the case of excitation by a single pulse, 
which turns out not to produce squeezing, and then turn 
to a two-pulse excitation, where under suitable condi- 
tions squeezing is found. Finally, in Sec. IIVI we finish 
with some concluding remarks. 



II. MODEL SYSTEM 

We consider a QD in the strong confinement limit. As- 
suming excitations by circularly polarized light, we can 
restrict the electronic states to a two-level system consist- 
ing of the ground state \g) and the lowest exciton state 
\x). These states are coupled to a classical light field as 
well as to LA phonons. The Hamiltonian of the system 
then reads 

H = h[n + 52(g^ + g&)]\x)(x\ 
q 

+fi£a; q &t& q -P-E, (1) 
q 

where HQ, is the exciton energy, 6 q (b q ) denotes the cre- 
ation (annihilation) operator of a LA phonon with wave 
vector q, cj q = cq is the phonon dispersion relation with 
the longitudinal sound velocity c, and g q is the electron- 
phonon coupling matrix element. We restrict ourselves 
to the case of deformation potential interaction, which 
for typical InAs/GaAs quantum dots has been found to 
be the dominant interaction mechanism on a picosecond 
time scaled Assuming for simplicity a spherical QD ge- 
ometry and a harmonic oscillator confinement, the cou- 
pling matrix element reads 
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with g being the crystal density and V the normaliza- 
tion volume of the crystal. D e (D h ) are the defor- 
mation potentials of electrons (holes) and a e (ah) the 
spatial widths of the electron (hole) wave functions. 
The classical laser field E is coupled to the polarization 



P = M \x)(g\ + M.Q\g)(x\ with the dipole matrix ele- 
ment Mo. We consider ultrafast laser pulses that are, 
however, spectrally narrow enough to realize a selective 
resonant excitation of the exciton. Once the selectivity 
has been accounted for by keeping only the resonantly 
coupled electronic levels the pulse duration is the short- 
est time scale in the problem and we can safely model 
E(t) as a series of delta- functions. We take GaAs ma- 
terial parameters and a QD with L = a e 2\/ln 2 = 5nm 
diameter. 26 

A basic quantity for the lattice dynamics is the lattice 
displacement associated with LA phonons 



u(r) = -i 
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In this paper we are particularly interested in the fluc- 
tuation squares of the lattice displacement, (Au(r)) 2 = 
(u(r) 2 ) — (u(r)) 2 . Assuming a spherical QD, all quantities 
depend only on the distance r from the QD center. For 
the interpretation of the results we introduce the relative 
excitation induced fluctuations squares 



D u (r,t) = 



(Au(r,t)) 2 -(Au vac ) 2 
(Au vac ) 2 



(4) 



where (A^ vac ) 2 are the fluctuation squares of the lattice 
displacement in the phonon vacuum state. For simplic- 
ity we will refer to D u just as fluctuations in the follow- 
ing. We call a state "squeezed" if (Au) 2 < (Aii vac ) 2 . 
Thus squeezing manifests itself directly by fluctuations 
D u (r) < 0. For LA phonons the vacuum fluctuation 
(Au vac ) 2 is calculated by taking the definition accord- 
ing to Eq. (j3j) and integrating over the whole Brillouin 
zone assumed as spherical as in the Debey model, i.e., 
up to the maximal wave vector g m ax = (247r 2 ) 1 / 3 /a for 
a zincblende structure with the lattice constant a. One 
finds 



(A^ vac ) 2 
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5.83- 10" 6 nm 2 . (5) 



For completeness we mention that the vacuum uncer- 
tainty of the lattice momentum (A7r vac ) 2 for LA phonons 
is given by 

(A^c) 2 = ( ^) 4/3 « 2.09 • 10 4 ( ^ " 

16 \ 7T / V nm 

(6) 

showing that for LA phonons the uncertainty product is 
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(A^ vac ) 2 (A7r vac ) 2 = - — 
o 4 



(7) 



and thus exceeds the Heisenberg limit already in the 
ground state. 

To calculate the expectation values of the lattice dis- 
placement and its fluctuations, we need expectation val- 
ues of the types (b q ), (b q b q >) and (b q b q >). For the case 
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of excitation by ultrafast laser pulses analytical results 
for these quantities can be found within a generating 
function formalism^ which will form the basis for the 
analytical expressions given in Sec. IIII1 



III. RESULTS 



A. Single pulse excitation 




Figure 1: (Color online) (a) Expectation value (u) -r and (b) 



fluctuations D u 



as functions of time t and distance r from 



the QD center after the excitation with •< 
pulse of pulse area 7r/2 at t = 0. 



single ultrafast laser 



Let us start by looking at the phonon dynamics after 
excitation by a single ultrafast laser pulse. In our previ- 
ous study of the fluctuation dynamics of LO phonons 12 
we have seen that the most interesting phenomena ap- 
pear in the case of excitation by pulses with a pulse area 
of 7r/2. It turns out that the same holds for LA phonons. 
Therefore, in this paper we will restrict our analysis to 
such pulses. We want to mention, that analytical results 
for the phonon dynamics like the ones shown here can be 
derived for any pulse area. 

An excitation with a single laser pulse of pulse area 
7r/2 creates an equal superposition of ground state \g) 
and exciton \x) in the electronic part of the system. The 
creation of the exciton causes a change in the charge car- 
rier distributions in the QD. Due to the deformation po- 
tential interaction this gives rise to a local deformation of 
the lattice corresponding to a shift of the equilibrium po- 
sitions of the lattice ions. In the case of acoustic phonons, 
this instantaneous shift in the equilibrium positions leads 
to the creation of a localized polaron in the QD and a 
wave packet leaving the QD^ The expectation value of 
the lattice displacement after the pulse reads 

(u(r,t)) = uoU(r,t) = u [W(r,t) + P(r)] (8) 

with uq = 8 2 2 • U(r,t) denotes the normalized lattice 
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displacement created by a single pulse excitation which 
can be decomposed into the part W(r, t) describing the 



wave packet leaving the QD given by 
1 



W(r,t) 



r I a e 



fZZ (r-ct 



7TT 
1 

r l 
2 

7rr 2 



D e erf 



r — ct 



ah 



D h erf 



(9) 



r — ct 
ah 



HZ (r+ct) 2 

D eVl e —^- 

a e . 



/— (r+ct) 2 

D h —e 

ah 



D e erf( I^)_^erf /r + C * 



and P(r) = —W(r,t = 0) describing the localized po- 
laron in the QD. Here, erf(x) denotes the error function. 
For the wave packet the dominant part for t ^> L/c is the 
first line in Eq. ([9]) that is a Gaussian centered around 
r = ct and decaying ~ 1/r. The other terms are decay- 
ing with 1/r 2 or describe an incoming wave packet that 
only contributes for small t. For small t, however, and 
thus also for the calculation of the polaron part P(r) all 
terms of Eq. (J9j) must be taken into account, because the 
divergences in front of the Gaussians and the error func- 
tions compensate each other resulting in a finite value at 
r = 0. 

Figure H{a) shows the expectation value of the lattice 
displacement (u) multiplied by r plotted as a function of 
time t after the laser pulse and distance r from the QD 
center. In this Figure the polaron can be identified as 
the horizontal line at r w 4 nm. The phonon wave packet 
that leaves the QD with the sound velocity c w 5nm/ps 
is seen as diagonal line. The amplitude of both polaron 
and wave packet are negative. 

In Fig. QJb) the corresponding fluctuations D u multi- 
plied by r 2 are plotted as a function of time t and po- 
sition r. It is clearly seen that D u is restricted to the 
same space-time regions where also the displacement is 
non- vanishing. Thus, the two lines in the plot can be 
identified as fluctuations of the polaron and of the trav- 
eling wave packet. Respectively D u has the same shape 
as the displacement. This is confirmed by the calcula- 
tions showing that for a single pulse excitation there is 
always D u ~ (u) 2 . Because D u is positive all the time 
we find that, as in the case of LO phonons, a single ultra- 
fast, resonant excitation of the QD never creates squeezed 
phonons in this system. 



B. Two pulse excitation 

In the next step we consider an excitation with two 
laser pulses each with a pulse area of tt/2. The first 
pulse arrives at time t = — r (with r > 0). The second 
pulse arrives at t = and has a relative phase of <p with 
respect to the first pulse. The expectation value of the 
displacement for t > after this two-pulse excitation can 
be written as 

(fl(r, t))=u [U(r, t + r) + Re (B~ (r)e^) U{r, t)] , (10) 
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Figure 2: (Color online) Real and imaginary part of the in- 
teraction amplitude B" as a function of the delay time r. 



where U (r, t) is the normalized lattice displacement after 
a single pulse [see Eqs. (|5J), ©]• The first pulse excit- 
ing the system at t = — r gives rise to the displacement 
U(r, t+r), the second pulse arriving at t = then creates 
U (r, £). The two terms are connected via the interference 
amplitude 



B (r) = exp 



(11) 



and the phase factor e^. For the coupling of Eq. ([2]) 
the sum in Eq. ([TT]) may be performed analytically (cf. 
appendix) . 

The real and imaginary part of B~ as functions of r 
are shown in Fig. [2j The real part exhibits an initial 
decay on the time scale of about 1 ps and then saturates 
at a value slightly above 0.9. The imaginary part is only 
nonzero during the initial decay of the real part. The 
interference amplitude is the same function that describes 
the decay of the optical polarization induced by the first 
laser pulse. 17 This decay is caused by the traveling wave 
packet which after about 1 ps has left the QD resulting in 
an entanglement of the QD with its environment and thus 
to a decoherence in the QD degrees of freedom. Since 
the second pulse couples to the polarization of the QD 
resulting from the first pulse, this decoherence reduces 
the coupling and therefore also reduces the generation of 
phonons by this second pulse. 

We want to mention that the saturation value of 
roughly 0.9 is also in good agreement with experimen- 
tal results obtained from time-integrated four- wave mix- 
ing signals on ensembles of QDs, 25 where an initial decay 
of the signals of about 20% at a temperature of 4 K has 
been found. According to the theory of four- wave mixing 
signals for the present type of models 2 ^ this initial decay 
is given by 1 — \B~(oo)\ 2 . This confirms the choice of the 
parameters used in our present calculations as realistic 
ones. 

Figure [3] (a) shows the lattice displacement created af- 
ter excitation with a pair of tt/2 pulses arriving at times 
t = — 3ps and t = with a relative phase of <p = 0. 
We clearly see two emitted wave packets, one starting at 
t = — 3ps the other at t = 0. At the time of the second 
excitation the polaron amplitude is essentially doubled. 
For a better interpretation of the results we note that for 
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Figure 3: (Color online) Same as Fig. \T\ but after the excita- 
tion with two laser pulses with pulse areas ty/2 at t = — 3ps 
and t = with a relative phase of (f> = 0. 



delay times longer than 2ps the interference amplitude 
B~ is essentially real and its value is larger than 0.9. 
Therefore the amplitude of the displacement created by 
the second pulse is not much smaller than the one created 
by the first pulse. To simplify the discussion we set in the 
following B~ w 1, which corresponds to neglecting the 
dephasing of the polarization after the first pulse. Note, 
however, that all the results shown in the Figures have 
been calculated with the correct value of B~ . 
Using B~ « 1, Eq. ([TO]) becomes 

(u(r,t)) ~u o [U(r,t + T)+cos(0)U(r,t)]. (12) 

In the case of a phase difference <fi = as taken in Fig. [3] 
the expectation value then evaluates to 

(u(r, t)> 1 0=o ~ u [2P(r) + W(r, t + r) + W(r, *)] , (13) 

where P(r) and W(r,t) are defined in Eqs. (|8J) and ([9]). 
This can be understood, when we look at the dynamics 
of the electronic system. The second pulse excites the 
system from the equal superposition of \g) and \x) to the 
exciton state \x). This doubling of the exciton occupation 
causes a doubling of the polaron amplitude in the QD. 
Because this process is similar to that induced by the 
first pulse, an identical wave packet is emitted. 

The fluctuations after a two pulse excitation with a 
large delay assuming B~ w 1 read 



D u (r,t) 



y[(W(r, t + r) 



(14) 



U(r, t)(sin 2 (0)/7(r, t) + 2 sin(0)C(r, t; r))] . 



In the fluctuations interferences between the single exci- 
tation processes take place. They are summarized in the 
term C(r, t; r). Its detailed form is given in the appendix. 

In Fig. [3]^b) we have plotted the corresponding fluctu- 
ations D u • r 2 . Positive values are seen at the positions of 
the two wave packets. The horizontal line of the polaron, 
however, is almost missing after the second pulse. These 
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two features can directly be seen from the fluctuations 
which in this case read 



D u (r : t)\ 



0=0 



^[W(r,t + T) 2 + W(r,t) 2 }. (15) 



Here we also assumed that W(r, t + r) • W(r, t) ~ 0, be- 
cause the overlap of the two wave packets is negligible. 
The physical reason for the absence of additional fluctu- 
ations in the QD beyond the vacuum fluctuations is the 
fact that for a QD being completely in the exciton state 
the polaron corresponds to a multimode shifted vacuum 
which has the same fluctuation properties as the phonon 
ground state. So, in the present case the polaron has fluc- 
tuations D u of almost zero and the positive fluctuations 
corresponding to the two emitted wave packets remain. 
No squeezed states are thus created by this laser pulse 
sequence. 




Figure 4: (Color online) (a), (b) same as Fig. [3] but for a 
relative phase (j) = 3tt/2. (c) energy density pE(r,t) • r 2 . 



From Eq. ([14]) one can see that squeezing can only 
occur if sin(0)?7(r, t)C(r, t; r) < 0. For <\> = utt this term 
vanishes. Thus, if squeezing occurs at all, we expect it 
to happen most prominently at phase differences around 
odd multiples of tt/2. In particular we will now study 
the case (j) = 3n/2 where sin(0) = — 1. Calculating the 
expectation value of the displacement for this case, the 
outcome of Eq. ([T2]) is 



</>=3tt/2 



u U(r,t + r), 



(16) 



because cos(37r/2) = 0. In other words, the second ex- 
citation of the system does not change the mean lattice 
displacement at all. Again looking at the electronic sys- 
tem, the second pulse does not act on the occupation in 
the electronic part of the system, we note that it just 



changes the relative phase in the equal superposition of 
\g) and \x). This case is shown in Fig. Ufa), and indeed 
the expectation value looks exactly like that of a single 
pulse excitation in Fig. Ufa). 

It is interesting to note that according to the exact 
result for the displacement given in Eq. ([TO]) , even in 
the case of arbitrarily strong dephasing there is always a 
value of the relative phase given by 



tan(</>) 



Re(B-) 
Im(B-) 



(17) 



where the second pulse does not modify the mean dis- 
placement. 

In contrast, the fluctuations are strongly affected by 
the second pulse, as can be seen in Fig. H(b). Du • r 2 
clearly shows fluctuations where the second wave packet 
would be expected. Thus our analysis shows that there 
is an emission of a phonon wave packet also in this case, 
however, with a vanishing mean displacement. We ob- 
serve that in the leading part of this wave packet clear 
negative parts build up in D u . So a squeezed phonon 
state is created under these excitation conditions. 

Due to the non-zero fluctuations energy transport 
should occur even where the expectation value (u) is zero. 
This energy can in principle be detected in bolometric 
measurements ] 23 ! 24 To quantify this, we calculate the en- 
ergy density of the LA phonons given by^ 



PE(r,t) = ^XV' 



2V 



) q UJ q > 



q,q 



(18) 



x Re 



i(q-q')-r 



(S q 6 q /)e 



i(q+q')-r 



where the vacuum energy density has already been sub- 
stracted. Again, after integrating over the Brillouin 
zone, the energy density only depends on the distance 
r from the QD. We present p#(r, t) • r 2 in Fig. H^c). 
Corresponding to the emission of the first wave packet 
the energy density is non-zero around the diagonal line 
(—3,0) — >• (2,30). Also where the second wave packet 
is expected, namely at the diagonal line (0, 0) — >• (6, 30) 
the energy density is non-zero. Here it is mainly positive, 
and a small negative value appears, corresponding to an 
energy below the vacuum energy. 

So far we have considered excitations, where the emit- 
ted wave packets are well separated. When we reduce 
the delay to r = 0.5 ps, the wave packets created by the 
two laser pulses are overlapping in space. We take again 
a relative phase of (j) = 3tt/2 which resulted in the ap- 
pearance of squeezing in the emitted wave packet. For 
this case D u is shown in Fig. [5fb). The mean displace- 
ment field is shown in Fig. G3a), it again does not exhibit 
much differences compared to the single pulse case. This 
can be understood from the fact that also here Im(B~) 
is much less than Re(B~) and thus Eq. ([T7]) is still well 
satisfied for <\> = 3tt/2. 

In the fluctuations D u we now find a broad wave packet 
with positive values in the leading and the trailing part 
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Figure 5: (Color online) Same as in Fig. [4j but with the first 
excitation at t — — 0.5ps. 



and pronounced negative values in between. The interac- 
tion of the phonons created by the second pulse with the 
wave packet resulting from the first pulse, which has not 
yet left the QD region, enhances the squeezing almost by 
a factor of two compared to the case of separated wave 
packets. 

The reduction of the delay time has a similar effect on 
the energy density, as is shown in Fig. [5fc). Here two 
significantly negative parts build up between the positive 
parts at the leading and trailing edges as well as at the 
center of the wave packets. 
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Figure 6: (Color online) Minimum values of D u -r 2 at t = 6ps 
and 20 nm < r < 50 nm as function of phase (j) and delay r. 

In order to obtain a complete picture of the squeezing 
behavior, when varying either the delay time r or the rel- 
ative phase (j) between the two laser pulses, Fig. [6] shows 
the minimum value of D u • r 2 in the emitted wave packets 
plotted as functions of r and (j). We find that only rela- 
tive phases in the range tt < <p < 2tt lead to remarkable 
squeezing. For < <fi < tt the achievable squeezing values 



are much smaller. When looking at the r-dependence of 
the plot, the largest fluctuations D u are obtained for de- 
lay times between 0.5 ps and 1 ps. In this range the wave 
packets generated by the two laser pulses exhibit a spatial 
overlap which obviously favors the build up of squeezing. 
For delays r w 1 ps the two wave packets have a spatial 
distance of r = cr ~ 5 nm, which is approximately the 
size of a single wave packet. For larger delays the spatial 
overlap is small and the squeezing values are reduced. 



IV. CONCLUSIONS 

In conclusion, we have analyzed the fluctuation prop- 
erties of LA phonons after ultrafast excitation of a QD 
in resonance with the lowest exciton transition. Like in 
the case of LO phonons ^ 11 ! 12 also for LA phonons we have 
found that a single pulse never gives rise to squeezing, but 
a sequence of two such excitations can create squeezed 
phonons. In contrast to the optical phonons, which due 
to their vanishing group velocity remain confined in the 
QD, the excitation of LA phonons leads to the formation 
of wave packets that leave the QD. We have shown that 
the phononic wave packet emitted after the excitation 
by the second pulse can be squeezed. In contrast to the 
squeezing found for LA phonons generated from the de- 
cay of LO phonons^ here also the phonon wave packets 
that leave the dot exhibit squeezing. The appearance of 
squeezing depends crucially on the relative phase (j) be- 
tween the two pulses. Sizeable squeezing is only found 
in the range ir < (j) < 2tt. For short delay times, when 
there is an overlap between the wave packets generated 
by the two pulses, the interaction between these wave 
packets has a significant impact on the strength of the 
fluctuations. 

Let us finally comment on the achieved values of the 
squeezing, which at first seems to be rather small. How- 
ever, it should be noted that vacuum fluctuations at a 
given point result from all phonon modes in the first Bril- 
louin zone. Excitation induced fluctuations, on the other 
hand, occur only for those phonon modes which couple 
to the QD exciton. According to the coupling matrix el- 
ement of Eq. ([2]) these are only phonons with vectors up 
to about the inverse size of the QD. If we would consider 
in the calculation of D u only those phonon modes which 
couple to the QD exciton a rather pronounced squeezing 
would emerge. Such a reduction of the effective range of 
phonon modes occurs also in an experiment if the mea- 
surement of the fluctuations is performed with a finite 
spatial resolution, which again eliminates the fluctuations 
resulting from phonons with wave vectors larger than the 
inverse spatial measurement resolution, 14 



Appendix A 

For completeness here we show the full analytical forms 
of some quantities introduced in the main text. First the 
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fluctuations for the two pulse excitation from Eq. (fT4|) 
without approximation read 



Du(r,t) = -JL[( W (r,t + T)-W(r,t)) 2 (Al) 



U(r, t)((l - R 2 )U(r, t) + 2IC(r, t; r))] 



is a generalized hypergeometric function. 

The complete analytical form of the energy density af- 
ter integration over q and q' is 



with R = Re(B-e^) and I = lm(B-e i(f) ). Here and in 
Eq. ([TTj) the interference amplitude appears; its full form 
is: 



B = exp 



exp 



V r 

2^ Jo 



l)dq 



(A2) 



47T 3 / 2 ghc 2 ' 



(A3) 



D e2 

2aJ 



i + erfi ( — ) 

D e D h 

>H^) 3/2 



2^ L 



z + erfi ( — ) 
Q>hJ 



i + erfi 



cr 



PE(r,t) 
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+ £(r, - 1)5 (r, t + r) - J£ (r, t; r))] 

£(r,t) = Ex(r - ct) + Ex(r + ct) - 2E 1 (r) 

+ E 2 {r-ct)-E 2 {r + ct), 

£(r,t) = 5i(r,t;r)+5 2 (r,t;r), 

5i(r,t;r) = #i (r - ct) + £?i (r + ct) 

- Ii(r-c(t + r))-li(r + c(t + r)), 
£ 2 (r,t;r) = E 2 (r - ct) + E 2 (r + ct) 

- E 2 (r - c(t + r)) - £ 2 (r + c(t + r)) 



(A4) 



with the complex error function erfi(x) = i ■ erf(ix). The 
interference term itself from Eqs. (|A1[) and (fl4|) is: 



C{r,t-T) 



with 



C(r - ct) 



and 



C(r - ct) + C(r + ct) 

[C(r-c(t + r)) + C(r + c(t + r))] 



/7T 



ct\ _(r-cty 
e a e 



e a l 



1 

2^ 



at 



(r - ct) 5 



H 



(r - ct) 



2 \ 1 



(r - ct) 2 



E^r-ct) 
E 2 (r-ct) 



-e a e 



E 1 (r-ct) 



— — e a e e a n 

t 1 \ a e a h 



e a h (r — ct) 



D h (r-cb\ -*=$£\ ( 
— erfi ( ) e a h (r — ct) 



E 2 (r-ct) = Ei(r-ct) 



\ a, 



■ erfi 



_(r-ct)^ 

e a £ 



#(z)= 2 F 2 (l,l; -,2;x 



oo 



2x™ 



S (2„+i)!!.( n +l) 
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